Transition from molecular complex to quantum solvation in ^HcivOCS 
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We present quantum calculations of the rotational energy levels and spectroscopic rotational 
constants of the linear OCS molecule in variable size clusters of ''He, using spectral evolution quan- 
tum Monte Carlo methods that allow excited states to be accessed without nodal constraints. The 
rotational constants of OCS are found to decrease monotonically from the gas phase value as the 
number of helium atoms increases to = 6, after which the average constant increases to satura- 
tion at a value in excellent agreement with experimental measurements made on significantly larger 
clusters (A'^ > 1000). The minimum is shown to indicate a transition from a molecular complex to a 
quantum solvated molecule, with the former characterized by floppy but near rigid behavior, while 
the latter is characterized by non-zero permutation exchanges and a smaller extent of rigid coupling. 

PACS Numbers: 36.40.-c, 05.30.Jp, 61.46. 4-w, 67.40.Yv 



Helium droplets offer a unique opportunity to immerse 
molecules into a superfluid. Spectroscopic experiments 
on molecules embedded in these droplets have yielded 
an array of measurements that provide on the one hand 
access to elementary excitations of the finite quantum 
liquid, and on the other hand information about the 
novel solvation dynamics encountered in a superfluid [Q . 
The rotational dynamics of the linear OCS molecule has 
played a key role in these spectroscopic studies. When 
solvated by the bosonic ^He isotope at temperatures 
T < 0.5 K, infra-red spectra of OCS show rotational fine 
structure that is consistent with a free rotation accompa- 
nied by an increased molecular moment of inertia, while 
no fine structure is observed when OCS is solvated by 
fermionic '^He clusters Quantum Monte Carlo (QMC) 
calculations have confirmed that a spectrum of free rota- 
tional states with modified moment of inertia will result 
from bosonic solvation Such free rotation is not just 
the result of solvation by a fluid interacting with weak 
van der Waals forces. It has been shown to be consistent 
with a negligible transfer of angular momentum from the 
molecule to relative motion of the solvating helium when 
the latter possesses Bose permutation symmetry . In 
contrast, efficient coupling to single-particle excitations 
in the fermionic '^He droplets has been shown to consider- 
ably reduce the lifetime of the rotational excitations [Q , 
accounting for the absence of rotational spectral transi- 
tions. 

Calculations of excited states afforded by QMC yield 
not only a direct route to calculation of spectroscopic 
constants for molecules solvated in helium, but they also 
provide critical insight into the nature of these rotational 
excitations. They thereby allow detailed microscopic 



analysis of the quantum coupling between a molecule 
and its helium environment. Direct calculations of ro- 
tational energy levels have been made within dynamical 
approximations such as the quasi-adiabatic rigid coupling 
approximation |^,|| and fixed node approximations [^,0 . 
However, in order to provide a complete understanding 
of the underlying quantum dynamics, it is necessary to 
use a theoretical approach that allows exact calculation of 
excited states, devoid of any dynamical or nodal approx- 
imations. This is possible with the projection operator 
imaginary time spectral evolution (POITSE) approach, 
a diffusion-Monte-Carlo-based methodology that has re- 
cently been applied to calculation of rotational excita- 
tions of the linear HCN molecule in helium clusters 0. 
For HCN, a considerably lighter molecule than OCS, 
these results showed a complex coupling between molec- 
ular and helium degrees of freedom that was attributed 
to the high zero point energy of the He-HCN system. In 
contrast, for the heavier octahedral molecule SFg, rota- 
tional constants derived from fixed node calculations em- 
ploying the free molecule nodal surfaces showed excellent 
agreement with experimental measurements 

In this letter we investigate the rotational motion of 
the OCS molecule inside ^He clusters, with direct calcu- 
lations of rotational energy levels by POITSE. Analysis of 
the results in terms of rigid coupling approximations and 
of permutation exchanges visible in related finite temper- 
ature calculations allow for the first time a clear transi- 
tion from a molecular complex to a molecule solvated by 
a quantum liquid to be identified as a function of clus- 
ter size N. OCS is intermediate in mass between HCN 
and SFg. Comparison of the rotational behavior of OCS 
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with that of HCN and of SFg provides key insights into 
the quantum couphng of the molecular rotational mo- 
tion with the solvating helium degrees of freedom, as a 
function of both molecular mass and symmetry. 

Our calculations employ the Projector Operator Imag- 
inary Time Spectral Evolution (POITSE) method for 
calculation of excited states In this scheme, ex- 

cited state energies are extracted from the two-sided 
inverse Laplace transform of an imaginary time corre- 
lation function k(t) that is computed by combining a 
multi-dimensional Monte Carlo integration with Diffu- 
sion Monte Carlo sidewalks. The decay of the correla- 
tion k(t) contains information about energy differences 
Ef — Eq, where Eq is the ground state energy and Ef 
an excited state energy level. The POITSE approach is 
ideal for calculation of excitations in many-particle sys- 
tems when estimates of excitation functions that take the 
ground state to the desired excited states are available. 
A suitable correlation function is then 

{^T \ Aexp[-{H - Eo)t]A^ \ ^t) ... 

k(t) — , (il 

{^T I exph(i7-So)T] I *t) 

where A is a local operator that projects from a trial func- 
tion I '^'t) approximating the ground state | ^to). The 
resulting initial state A' \ '^'t) is usually not an eigen- 
function of H, but time evolution under the action of the 
full Hamiltonian will ensure that all cigcn components 
4*/) that overlap with this state will contribute to the 
correlation function. Inverse Laplace transform of k(t), 
performed by the maximum entropy method as described 
in Ref. 1^, results in the spectral function 

/ 

from which all excitation energies accessed by A may be 
extracted as the peak locations. When | ^>t) =| 5'o)> 
exact energies are obtained, while otherwise some trial 
function bias may exist 

Our approach to extract the rotational excitations of 
interest is to take a free molecular projector acting on 
the molecule-helium cluster ground state. Since the 
molecule-helium coupling is weak, this provides a good 
starting point for calculation of the imaginary time corre- 
lation function and subsequent spectral transformation. 
For the trial function ^>t we employ variationally op- 
timized ground state functions of the usual generalized 
product form, namely containing pairwise correlations 
between all components of the cluster, i.e., isotropic He- 
He and anisotropic He-OCS terms. The He-He correla- 
tions are those used in Ref. while the He-OCS corre- 
lations take the form x{R, {pi}) = exp{po-R^^ +P2[1 + 
P3 cos{9 — Pi) \nR] + p^R^ cos^(6' — p4)e(P<'~^'^^)} where 
_R, 6 are the Jacobi coordinates describing the location 
of the He atom relative to the linear OCS molecule [p| . 



We consider here projectors A proportional to the molec- 
ular Wigner functions i?^^. in a space-fixed frame. Ze- 
ros of Dl^f^ constitute free molecule nodal surfaces for 
a symmetric top [ pl]| . We focus here in particular on 
j = l,m = fc = 0,a linear rotor eigenfunction for which 
A=cos /3 is a function only of the second Eulcr angle of 
the molecule specifying the orientation of the molecular 
frame in the (arbitrary) space-fixed frame. A accesses 
states in which the total angular momentum J is carried 
primarily by the OCS, i.e., the molecular angular mo- 
mentum j is a quasi-good quantum number and there is 
negligible angular momentum / in the relative helium mo- 
tion (J = j -I- 1). Our choice of Dqq is motivated by analy- 
sis of the calculated rotational spectrum for the He-OCS 
complex {N — 1) which shows that this is the pri- 
mary molecular function contributing to the JKpKa=^m 
rotational state of the complex. (We employ here the 
notation J, Kp, Kq for rotational states of an asymmetric 
rotor 0). 

We emphasize that the projector A provides only an 
initial guess for the nodal structure of the excitations, 
and that these automatically adjust to their true val- 
ues during the calculation. That these arc indeed sub- 
tly different from the free nodal surfaces for OCS is evi- 
dent from the fact that for = 1, a 5-dimensional sys- 
tem where exact calculations are possible using e.g., the 
BOUND program jl^ , fixed node calculations carried out 
with nodal constraints imposed by the Wigner function 
do not show good agreement with the exact calculations, 
overestimating the exact rotational excitation energy for 
OCS-He (0.307 cm-^) by 0.07±0.01 cm-^, i.e., 22(3)%. 
In contrast, POITSE gives the exact energy to within a 
statistical error of 0.7%, i.e., 0.309±0.002 cm'^. This 
behavior contrasts with the high accuracy of the free 
molecule fixed node approximation found for the heav- 
ier SFg-He complex. ^ 

The total molecule-cluster interaction potential is the 
sum of all pairwise contributions (He-He and OCS- He). 
We have employed two recent He-OCS potentials [ p^jl3 
and find equivalent results with these. Unless otherwise 
stated, all calculations reported here are made with the 
ah initio HHDSD He-OCS potential energy surface |l^] 
and the HFD-B He-He potential of Ref. The imag- 
inary time evolution is performed with the rotational 
importance-sampled rigid body Diffusion Monte Carlo 
algorithm of Ref. Initial ensembles of 1000 walk- 

ers distributed according to \^t\^ were propagated for 
30,000 steps of At=50 a.u. using the mixed branch- 
ing/weight algorithm of Ref. |l§. Typically 1000-2000 
independent decays were required to produce a converged 
spectrum k{u)), with the largest size (TV = 20) requiring 
the most decays. Statistical error bars in the excitation 
energies were estimated using a Gaussian approximation, 
according to the procedure described in Ref. . 

POITSE calculations with the Wigner projector A 
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consistently produced only one peak in the k{uj) spec- 
trum, for all sizes N < 20. This clean single-valued 
behavior of the free molecule projector for OCS con- 
trasts with the behavior seen earlier for the lighter HCN 
molecule, where the same projector leads to multiple ex- 
cited states Q. Single- valued free molecule projectors 
suggest that AI'^/t) is very similar to the true eigenstate. 
However, free or quasi-free rotations do not prohibit ef- 
fective rigid coupling of some solvating helium density to 
the molecular rotation, leading to a renormalization of 
the molecular moment of inertia similar to that found in 
the microscopic two-fluid theory |^ and in perturbative 
analysis 

When interpreted with an asymmetric top Hamilto- 
nian the cluster rotational excited state energy de- 
rived from A, E(Jk ,Ka = loi)i yields an effective rota- 
tional constant Bavg- When close to the limit of a sym- 
metric top, Bavg is equal to either [B+C)/2 (prolate top) 
or {A-\-B)/2 (olDlate top). Fig. |l| shows the main result of 
this paper, namely the behavior of Bavg obtained from 
the POITSE calculations, as a function of cluster size 
N . Corresponding experimental values of B for the free 
(gas phase) OCS molecule and for OCS in large helium 
droplets {N > 1000) where the spectrum was fit to a lin- 
ear rotor {Bavg = B) are shown from Ref. For the 
smallest clusters iV = 1 Q and iV = 2 we show the 
experimentally measured values of Bavg ■ 

The POITSE results show a monotonic decrease of the 
excited state energy up to = 6, followed by a small con- 
tinuous increase to saturation at a value Bavg = 0.07(2) 
cm~^ that is in excellent agreement with the experi- 
mentally measured value B = 0.0732(3) cm~^ in large 
droplets Excellent agreement of Bavg with the cor- 
responding experimental values Bavg for small clusters 
is also achieved for the two sizes A^=l, 2 for which this 
experimental data is available [|l3|,0. Thus the POITSE 
calculations are able to bridge the gap between small van 
der Waals clusters and large droplets, achieving high ac- 
curacy in both regimes. The overall behavior of Bavg for 
OCS shows some similarities with SFg and with HCN, 
but also exhibits critical differences from these molecules. 
For the more symmetrical and heavier SFg molecule the 
saturation value is also reached before the first solva- 
tion shell is completed {N ~ 8) but negligible overshoot 
is seen Q. For the lighter HCN, also linear, a large 
overshoot is seen but the rotational constant does not 
reach saturation before completion of the first solvation 
shell [|. 

A deeper understanding of the size dependent quan- 
tum coupling within these rotational excitations in the 
bosonic '^He cluster is provided by analysis in terms of 
both rigid coupling models and permutation exchange 
distributions. The triangle symbols in Fig. ^ show 
the values of Bavg obtained using the quasi-adiabatic 
rigid coupling approach ||^,^ implemented with corre- 
lated sampling j20|. The POITSE results are very close 



to the rigid coupling estimates for iV < 6, with the latter 
lying consistently slightly lower than the corresponding 
POITSE values. Within a rigid coupling approximation, 
the complete ^He density is assumed to adiabatically fol- 
low the rotational motion of OCS. Therefore the system- 
atically higher values of the POITSE results is a direct 
manifestation of only partial adiabatic following of the 
solvating helium density with OCS rotation. The quasi- 
adiabatic rigid coupling approach does incorporate all 
zero point motions. This results in the intersection at 
iV = 4 of the rigid coupling estimate of Bavg with the 
experimental B value in large droplets, in contrast to 
the crossing at A^=6 estimated in Ref. from a clas- 
sical binding model that neglects zero point energy. The 
actual crossing point is in remarkable agreement with 
the integral of the molecule-induced local non-superfluid 
density calculated by path integral methods in Ref. j|] 
{iT-ns ^ 3.2), and is thus consistent with the analysis 
made there of a reduced B deriving from rigid coupling 
to the local non-superfluid density in the first solvation 
shell. The quasi-adiabatic rigid coupling estimates fail at 
larger sizes, where they necessarily continue to decrease 
monotonically, reaching negligible values by TV 50. Ta- 
ble I shows the quasi-adiabatic rigid coupling results for 
all three rotational constants A, B and C, for selected 
cluster sizes. This indicates that the asymmetric top 
spectrum {A> B > C) evolves into a prolate symmetric 
top spectrum [A > B = C) hy N = 10, the same size at 
which the POITSE Bavg reaches the asymptotic droplet 
value of B (Fig. |l]). However, we expect that the transi- 
tion between asymmetric and symmetric top may occur 
at a different cluster size for the exact calculations. 

One of the most interesting features of the POITSE 
results seen in Fig. |l| is the turn-around of Bavg at rela- 
tively small N , and subsequent rise to saturation at the 
experimental large droplet value by iV = 10. This be- 
havior is related to the onset and nature of permutation 
exchanges between the ^He atoms. It can be quantita- 
tively interpreted in terms of these exchanges pl| and 
of the structures previously analyzed in terms of the 
axial solvation locations in Ref. |2^. As shown there, 
while for < 5 the helium density is essentially lo- 
calized in the global minimum of the He-OCS potential 
energy surface, as N increases beyond 5 the additional 
helium atoms solvate other regions along the molecule, 
covering the entire axial extent of the OCS molecule by 

= 10. Path integral calculations show no permutation 
exchanges for < 5, but as A^ increases above 5 and the 
helium density grows along the molecular axis, exchange 
permutations are seen, with components both along and 
around the molecular axis ||2l[ . The lack of exchanges for 
Af < 5 is attributed to the difficulty of exchanging within 
a single tightly packed axial ring, and accounts for the 
near rigid behavior at these sizes. The permutation ex- 
change path components along the molecular axis seen 
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for N >6 result in a lowering of the corresponding rigid 
body response for rotation about axes perpendicular to 
the molecule (i.e., to a non-zero perpendicular superfluid 
response) , causing in a rise in rotational constants B and 
C . This effect continues as the solvation layer grows 
along the molecular axis and the lateral permutation ex- 
change contributions increase, thereby accounting for the 
observed increase in the POITSE values over the range 
iV = 6 - 10 in Fig. Note that = 10 is the first size 
at which the *He density extends along the entire molec- 
ular axis pO| ] , thereby allowing permutation exchanges of 
maximal extent within the first solvation shell, i.e., from 
one end of the molecule to the other. We conclude that 
for OCS the range iV = 6 — 10 with its increase of Bavg 
to saturation from an overshoot at iV = 6 denotes a tran- 
sition between a van der Waals (molecular) complex and 
a true quantum solvated molecule. The complex is 'near 
rigid', or 'floppy', in spectroscopic language, while the 
quantum solvated molecule is characterized by permuta- 
tion exchanges between the helium atoms [||. For large 
enough clusters, the latter leads to a microscopic local 
two-fluid description of co-existing superfluid and non- 
superfluid densities within the first solvation shell p3|,pt. 

In conclusion, we have provided calculation of rota- 
tional constants for OCS-doped ''He clusters containing 
up to one complete solvation shell (N = 1 — 20), by a 
direct method without imposing dynamical or nodal ap- 
proximations. We find accurate agreement with experi- 
mental results at both extremes of small size N — 1,2, 
and large droplets. These POITSE calculations allow for 
the first time the transition from a molecular complex to 
a quantum solvated molecule in superfluid ^He to be pre- 
cisely identified, and interpreted in terms of the quantum 
structure and permutation exchange propensity of the lo- 
cal helium environment. 
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TABLE I. Rotational constants (in cm~^) calculated by 
the quasi-adiabatic rigid coupling method with two different 
He-OCS potential energy surfaces, shown for clusters having 
iV < 10 ''He atoms. A,B,C: from HHDSD potential @. 
A',B',C': from MP4 potential 



N 


A 


A' 


B 


B' 


C 


C 


1 


0.457(2) 


0.443(2) 


0.175(2) 


0.173(2) 


0.123(2) 


0.120(2) 


2 


0.208(2) 


0.202(2) 


0.128(2) 


0.129(2) 


0.092(1) 


0.092(1) 


4 


0.101(2) 


0.098(1) 


0.074(1) 


0.077(1) 


0.062(1) 


0.065(1) 


6 


0.067(1) 


0.065(1) 


0.051(1) 


0.052(1) 


0.044(1) 


0.045(1) 


7 


0.057(1) 


0.056(1) 


0.045(1) 


0.045(1) 


0.038(1) 


0.039(1) 
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FIG. 1. Effective rotational constant Bavg for OCS in 
*Hejv, as a function of cluster size A''. Circles show POITSE 
results obtained from cluster rotational state energies gener- 
ated by projectors A = cos(/3) (see text). Triangles show 
quasi-adiabatic rigid coupling estimates [M, asterisks corre- 
spond to the experimental values of Bavg for A^=l,2 
The latter are seen to be coincident with the POITSE results 
(circles). Gas phase and experimental values in large ^He 
droplets are taken from Ref. [hsj. 
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